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Weak gravitational lensing allows one to reconstruct the spatial distribution of the projected mass density 
across the sky. These “mass maps” provide a powerful tool for studying cosmology as they probe both luminous 
and dark matter. In this paper, we present a weak lensing mass map reconstructed from shear measurements in 
a 139 deg^ area from the Dark Energy Survey (DES) Science Verification data. We compare the distribution of 
mass with that of the foreground distribution of galaxies and clusters. The overdensities in the reconstructed map 
correlate well with the distribution of optically detected clusters. We demonstrate that candidate superclusters 
and voids along the line of sight can be identified, exploiting the tight scatter of the cluster photometric redshifts. 
We cross-correlate the mass map with a foreground magnitude-limited galaxy sample from the same data. Our 
measurement gives results consistent with mock catalogs from N-body simulations that include the primary 
sources of statistical uncertainties in the galaxy, lensing, and photo-z catalogs. The statistical significance of the 
cross-correlation is at the 6.8(7 level with 20 arcminute smoothing. We find that the contribution of systematics 
to the lensing mass maps is generally within measurement uncertainties. In this work, we analyze less than 3 % 
of the final area that will be mapped by the DES; the tools and analysis techniques developed in this paper can 
be applied to forthcoming larger datasets from the survey. 

PACS numbers: 


I. INTRODUCTION 

Weak gravitational lensing is a powerful tool for cosmo¬ 
logical studies [see 1, 2, for detailed reviews]. As light from 
distant galaxies passes through the mass distribution in the 
Universe, its trajectory gets perturbed, causing the apparent 
galaxy shapes to be distorted. Weak lensing statistically mea¬ 
sures this small distortion, or “shear”, for a large number of 
galaxies to infer the 3D matter distribution. This allows us to 
constrain cosmological parameters and study the distribution 
of mass in the Universe. 

Since its first discovery, the accuracy and statistical pre¬ 
cision of weak lensing measurements have improved signif¬ 
icantly [3-8]. Most of these previous studies constrain cos¬ 
mology through N-point statistics of the shear signal [e.g. 
9-14]. In this paper, however, we focus on generating 2D 
wide-field projected mass maps from the measured shear [15]. 
These mass maps are particularly useful for viewing the non- 
Gaussian distribution of dark matter in a different way than is 
possible with N-point statistics. 

Probing the dark matter distribution in the Universe is par¬ 
ticularly important for several reasons. Based on the peak 
statistics from a mass map it is possible to identify dark mat¬ 
ter halos and constrain cosmological parameters [e.g. 16-20]. 
Mass maps also allow us to study the connection between 
baryonic matter (both in stellar and gaseous forms) and dark 
matter [15]. This can be measured by cross correlating light 
maps and gas maps with weak lensing mass maps. Correla¬ 
tion with light maps, which can be constructed using observed 
galaxies, groups and clusters of galaxies etc., can be used 
to constrain galaxy bias, the mass-to-light ratio, and the de¬ 
pendence of these statistics on redshift and environment [21- 
24]. However, one needs to take caution when interpreting 
the weak lensing mass maps, as the completeness and purity 
of structure detection via these maps is not very high due to 
their noisy nature [25]. 

One other interesting application of the mass map is that it 
allows us to identify large scale structures (both super-clusters 
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and voids) which are otherwise difficult to find [e.g. 26]. Char¬ 
acterizing the statistics of large structures can be a sensitive 
probe of cosmological models. Structures with masses as high 
or higher than clusters require special attention as the massive 
end of the halo mass function is very sensitive to the cosmol¬ 
ogy [27-29]. These rare structures also allows us to constrain 
different theories of gravity [30, 31]. In addition to the study 
of the largest assemblies of mass, the study of number density 
of the largest voids allows further tests of the ACDM model 
[e.g. 32]. 

Similar mass mapping technique as used in this paper has 
been previously applied to the Canada-France-Hawaii Tele¬ 
scope Lensing Survey (CFHTLenS) as presented in Van Waer- 
beke et al. [33]. Their work demonstrated the potential sci¬ 
entific value of these wide-field lensing mass maps, includ¬ 
ing measuring high-order moments of the maps and cross¬ 
correlation with galaxy densities. The total area of the mass 
map in that work is similar to our dataset, though it was di¬ 
vided into four separate smaller fields. 

The main goal of this paper is to construct a weak lens¬ 
ing mass map from a contiguous 139 deg^ area in the Dark 
Energy Survey[89] [DES, 34, 35] Science Verification (SV) 
data, which overlaps with the South Pole Telescope survey 
(the SPT-E field). The SV data were recorded using the newly 
commissioned wide-field mosaic camera, the Dark Energy 
Camera [DECam; 36-38] on the 4m Blanco telescope at the 
Cerro Tololo Inter-American Observatory (CTIO) in Chile. 
We cross correlate this reconstructed mass map with optically 
identified structures such as galaxies and clusters of galaxies. 
This work opens up several directions for future explorations 
with these mass maps. 

This paper is organized as follows. In Sec. II we describe 
the theoretical foundation and methodology for constructing 
the mass maps and galaxy density maps used in this paper. 
We then describe in Sec. Ill the DES dataset used in this 
work, together with the simulation used to interpret our re¬ 
sults. In Sec. IV we present the reconstructed mass maps. We 
discuss qualitatively in Sec. V the correlation of these maps 
with known foreground structures found via independent op¬ 
tical techniques. In Sec. VI, we quantify the wide-field mass- 
to-light correlation on different spatial scales using the full 
field. We show that our results are consistent with expecta- 
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tions from simulations. In Sec. VII we estimate the level of 
contamination by systematics in our results from a wide range 
of sources. Finally, we conclude in Sec. VIII. For a summary 
of the main results from this work, see the companion paper 
in PRL [39]. 


II. METHODOLOGY 

In this section we first briefiy review the principles of weak 
lensing in Sec. II A. Then, we describe the adopted mass re¬ 
construction method in Sec. IIB. Finally in Sec. IIC, we de¬ 
scribe our method of generating galaxy density maps. The 
galaxy density maps are used as independent mass tracers in 
this work to help confirm the signal measured in the weak 
lensing mass maps. 


A. Weak gravitational lensing 

When light from galaxies passes through a foreground mass 
distribution, the resulting bending of light leads to the galaxy 
images being distorted [e.g. 1]. This phenomenon is called 
gravitational lensing. The local mapping between the source 
(/3) and image {6) plane coordinates (aside from an overall 
displacement) can be described by the lens equation: 

/3-/3o=A(0)(0-0o), (1) 

where (3q and Oq is the reference point in the source and the 
image plane. A is the Jacobian of this mapping, given by 

= (!-.)(■-*' -4), (2) 

where K is the convergence, gi = yi/{\ — k) is the reduced 
shear and Ji is the shear. i= 1,2 refers to the 2D coordinates 
in the plane. The factor (1 — /c) causes galaxy images to be 
dilated or reduced in size, while the terms in the matrix cause 
distortion in the image shapes. Under the Bom approxima¬ 
tion, which assumes that the deflection of the light rays due to 
the lensing effect is small, A is given by [e.g. 1] 

Aij{6,r) = (3) 

where if/ is the lensing deflection potential, or a weighted pro¬ 
jection of the gravitational potential along the line of sight. 
For a spatially fiat Universe, it is given by the line of sight 
integral of the 3D gravitational potential <I> [40], 

y/{e,r) = 2 (4) 

Jo rr ^ ^ 

where r is the comoving distance. Comparison of Eq. (3) with 
Eq. (2) gives 



(5) 


7 = ri+m 



¥,22) + i¥,l2- 


( 6 ) 


For the purpose of this paper, we use the Limber approxima¬ 
tion which lets us use the Poisson equation for the density 
fluctuation 5 = (A — A)/A (where A and A are the 3D density 
and mean density respectively): 

v2<D=^^^5, (7) 

2a 

where a is the cosmological scale factor. Eq. (4) and Eq. (5) 
give the convergence measured at a sky coordinate 0 from 
sources at comoving distance r: 


_ r 5 (6,/) 

2 Jo ^ r a{r') 


( 8 ) 


We can generalize to sources with a distribution in comoving 
distance (or redshift) f{r) as: k{6) = / K{6,r)f{r)dr. That 
is, a k: map constmcted over a region on the sky gives us the 
integrated mass density fiuctuation in the foreground of the K 
map weighted by the lensing weight p(r'), which is itself an 
integral over /(r): 


m = 





(9) 


with 


rrn r — r' 

p{r)= dr/(r)-, (10) 

Jr' r 

where vh is the comoving distance to the horizon. For a spec¬ 
ified cosmological model and /(r) specified by the redshift 
distribution of source galaxies, the above equations provide 
the basis for predicting the statistical properties of K. 


B. Mass maps from Kaiser-Squires reconstruction 

In this paper we perform weak lensing mass reconstruction 
based on the method developed in Kaiser and Squires [41]. 
The Kaiser-Squires (KS) method is known to work well up 
to a constant additive factor as long as the structures are in 
the linear regime [33]. In the non-linear regime (scales cor¬ 
responding to clusters or smaller structures) improved meth¬ 
ods have been developed to recover the mass distribution [e.g. 
42, 43]. In this paper we are interested in the mass distribution 
on large scales; we can therefore restrict ourselves to the KS 
method. The KS method works as follows. The Fourier trans¬ 
form of the observed shear, 7 , relates to the Fourier transform 
of the convergence, k through 

( 11 ) 


~ I /?I9 


( 12 ) 


where it are the Fourier counterparts for the angular coordi¬ 
nates 0 /, / = 1,2 represent the two dimensions of sky coor¬ 
dinate. The above equations hold true for i ^ 0. In practice 
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we apply a sinusoidal projection of sky with a reference point 
at RA=71.0 deg and then pixelize the observed shears with a 
pixel size of 5 arcmin before Fourier transforming. Given that 
we mainly focus on scales less than a degree in this paper, the 
errors due to the projection is small [33]. 

The inverse Fourier transform of Eq. (11) gives the con¬ 
vergence for the observed field in real space. Ideally, the 
imaginary part of the inverse Fourier transform will be zero 
as the convergence is a real quantity. However, noise, system- 
atics and masking causes the reconstruction to be imperfect, 
with non-zero imaginary convergence as we will quantify in 
Sec. VIB. The real and imaginary parts of the reconstructed 
convergence are referred to as the E- and B-mode of k, re¬ 
spectively. In our reconstruction procedure we set shears to 
zero in the masked regions [44]. We later quantify the effect 
of this step in Sec. VIB. 

One of the issues with the KS inversion is that the uncer¬ 
tainty in the reconstructed convergence is formally infinite for 
a discrete set of noisy shear estimates. This is because the sta¬ 
tistically uncorrelated ellipticities of galaxies result in a white 
noise power spectrum which integrates to infinity for large 
spatial frequencies. Therefore we need to remove the high 
frequency components. For a Gaussian filter of size a the co- 
variance of the statistical noise in the convergence map can be 
written as [45] 

where Gg is the standard deviation of the single component 
ellipticity (which contains the intrinsic shape noise and mea¬ 
surement noise) and rig is the number density of the source 
galaxies. Eq. (13) implies that the shape noise contribution to 
the convergence map reduces with increasing size of the Gaus¬ 
sian window and number density of the background source 
galaxies. 


C. Lensing-weighted galaxy density maps 

In addition to the mass map generated from weak lensing 
measurements in Sec. IIB, we also generate mass maps based 
on the assumption that galaxies are linearly biased tracers of 
mass in the foreground. In particular, we study two galaxy 
samples: the general field galaxies and the Luminous Red 
Galaxies (LRGs). Properties of the samples used in this work 
such as the redshift distribution, magnitude distribution etc. 
are described in Sec. IIIB. To compare with the weak lensing 
mass map, we assume that the bias is constant. However, bias 
may change with spatial scale, redshift, magnitude and other 
galaxy properties. This can introduce differences between the 
weak lensing mass map and foreground map. In this paper 
we neglect such effects since we mostly focus on large scales 
(>5 — 10 arcmin at z ^ 0.35) where the departures from linear 
bias are small [46]. 

Based on a given sample of mass tracer we generate a 
weighted foreground map (Kg) after applying an appropriate 
lensing weight to each galaxy before pixelation. In principle 


the weight increases the signal-to-noise (S/N) of the cross¬ 
correlation between the lensing mass map and the foreground 
density map. The lensing weight (Eq. (10)) depends on the 
comoving distance to the source and lens, and the distance 
between them. To generate the weighted galaxy density map, 
we first generate a 3D grid of the galaxies. We estimate the 
density contrast in each of these cells as follows: 


^ nk 


where (ij) is the pixel index in the projected 2D sky and k 
is the pixel index in the redshift direction, ntjk is the num¬ 
ber of galaxies in the i cell and fik is the average number 
of galaxies per pixel in the redshift bin. This 3D grid of 
galaxy density fiuctuations will be used to estimate Kg accord¬ 
ing to the discrete version of Eq. (9), 


jj _ . ^k^^k {dl dk)fl 


(15) 


where Kg is the weighted foreground map at the pixel (ij)', 
k and I represent indices along the redshift direction for lens 
and source, is the physical size of the redshift bin, J/ is 
the angular diameter distance to source, // is the probability 
density of the source redshift distribution at redshift I and 5^^ 
is the foreground density fiuctuation at angular diameter dis¬ 
tance dk. In this work, use a single source redshift bin and 
A^ = 0.1 for the lens sample. We adopt the following cosmo¬ 
logical parameters: = 0.3, = 0.7, Q.k = 0.0, h = 0.72. 

Our results depend very weakly on the exact values of these 
cosmological parameters. 


III. DATA AND SIMULATIONS 

The measurements in this paper are based on 139 deg^ of 
data in the SPT-E field, observed as part of the Science Veri¬ 
fication (SV) data from DES. The SV data were taken during 
the period of November 2012 - February 2013 before the of¬ 
ficial start of the science survey. The data were taken shortly 
after DECam commissioning and were used to test survey op¬ 
erations and assess data quality. Five optical filters (grizY) 
were used throughout the survey, with typical exposure times 
being 90 sec for griz and 45 sec for Y. The final median depth 
estimates of this data set in our region of interest are g ^ 24.0, 
r ^ 23.9, i ^ 23.0 and z ^ 22.3 (10-cr galaxy limiting magni¬ 
tude). 

Below we introduce in Sec. Ill A the relevant data used in 
this work. Then we define in Sec. IIIB two subsamples of the 
SV data that we identify as “foreground (lens)” and “back¬ 
ground (source)” galaxies for the main analysis of the paper. 
In Sec. me we introduce the simulations we use to interpret 
our measurements. 


A. The DES SVAl Gold galaxy catalogs 

All galaxies used for foreground catalogs and lensing mea¬ 
surements are drawn from the DES SVAl Gold Catalog 
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(Rykoff et al., in preparation) and several extensions to it. 
The main catalog is a product of the DES Data Management 
(DESDM) pipeline version “SVAl” (Yanny et ah, in prepa¬ 
ration). The DESDM pipeline, as described in Ngeow et al. 
[47], Sevilla et al. [48], Desai et al. [49], Mohr et al. [50], 
begins with initial image processing on single-exposure im¬ 
ages and astrometry measurements from the software pack¬ 
age SCAMP [51]. The single-exposure images were then 
stacked to produce co-add images using the software package 
SWARP [52]. Basic object detection, point-spread-function 
(PSP) modelling, star-galaxy classification [90] and photom¬ 
etry were done on the individual images as well as the co-add 
images using software packages SExtractor [53] and PS- 
PEx [54]. The full SVAl Gold dataset consists of 254.4 deg^ 
with gr/z-band coverage, and 223.6 deg^ for Y band. The main 
science goal for this work is to reconstruct wide-field mass 
maps; as a result, we use the largest continuous region in the 
SV data: 139 deg^ in the SPT-E field. 

The SVAl Gold Catalog is augmented by: a photometric 
redshift catalog, two galaxy shape catalogs, and an ERG cat¬ 
alog. These catalogs are described below. 


7. Photometric redshift catalog 

In this work we use the photometric redshift (photo-z) es¬ 
timated with the Bayesian Photometric Redshifts (BPZ) code 
[55, 56]. The photo-z’s are used to select the main foreground 
and background sample (see Sec. IIIB). The details and ca¬ 
pabilities of BPZ on early DES data were already presented 
in Sanchez et al. [57], where it showed good performance 
among template-based codes. The primary set of templates 
used contains the Coleman et al. [58] templates, two starburst 
templates from Kinney et al. [59] and two younger starburst 
simple stellar population templates from Bruzual and Char¬ 
iot [60], added to BPZ in Coe et al. [56]. We calibrate the 
Bayesian prior by fitting the empirical function n(z, ^ |mo) pro¬ 
posed in Benitez [55], using a spectroscopic sample matched 
to DES galaxies and weighted to mimic the photometric prop¬ 
erties of the DES SV sample used in this work. As tested in 
Sanchez et al. [57], the bias in the photo-z estimate is ^0.02, 
with 68% scatter 068 ^0.1 and the 3o outlier fraction ^2%. 
Eor this work, we use Zmean^ the mean of the Probability Dis¬ 
tribution Eunction (PDE) output from BPZ as a single-point 
estimate of the photo-z to separate our galaxies into the fore¬ 
ground and background samples. Other photo-z codes used in 
DES have been run on the same data. Eor this work we have 
also checked our main results in Sec. VI using an indepen¬ 
dent Neural Network code [Skynet; 61, 62]. We found that 
BPZ and Skynet gives consistent results (within 1 cr) in terms 
of the cross-correlation between the weak lensing mass maps 
and the foreground galaxy map. 


2. Shape catalogs 

Based on the SVAl data, two shear catalogs were produced 
and tested extensively in Jarvis et al. (in preparation): the 


ngmix [91] (version 011) catalog and the imSshape [92] (ver¬ 
sion 9) catalog. The main results shown in our paper are 
based on the ngmix catalog, but we also cross-check with the 
imSshape catalog to confirm that the results are statistically 
consistent. These catalogs are slightly earlier version from 
that described in Jarvis et al. (in preparation). 

The PSP model for both methods are based on the single¬ 
exposure PSP models from PSPEx. PSPEx models the PSP 
as a linear combination of small images sampled on an ad hoc 
pixel grid, with coefficients that are the terms of a second- 
order polynomial of pixel coordinates in each DECam CCD. 

ngmix [63] is a general tool for fitting morphological mod¬ 
els to images of astronomical objects. Eor the galaxy model, 
ngmix supports various options including exponential disk 
and de Vaucouleurs’ profile [64], all of which are imple¬ 
mented approximately as a sum of Gaussians [65]. Addition¬ 
ally, any number of Gaussians can be fit. These Gaussian fits 
can either be completely free or constrained to be co-centric 
and co-elliptical. Eor the DES SV galaxy images, we used 
the exponential disk model. Eor the PSP fitting, an Expecta¬ 
tion Maximization [66] approach is used to model the PSP as 
a sum of three free Gaussians. Shear estimation was carried 
out using by jointly fitting multiple images in r, /,z bands. The 
multi-band approach enabled a larger effective galaxy number 
density compared to the imSshape catalog, which is based on 
single-band images in the current version. 

The imSshape [67] implementation in this work estimates 
shapes by jointly fitting a parameterized galaxy model to all of 
the different single-exposure r-band images, finding the max¬ 
imum likelihood solution. Calibration for bias in the shear 
measurement associated with noise [68, 69] is applied. An 
earlier version of this code (run on the co-add images instead 
of single-exposures) has been run on the SV cluster fields for 
cluster lensing studies [70]. 

Details for both shape catalogs and the tests performed on 
these catalogs can be found in Jarvis et al. (in preparation). 
Both shear catalogs have been tested and shown to pass the re¬ 
quirements for S V cosmic shear measurement, which is much 
more stringent than what is required in this paper. As our anal¬ 
ysis is insensitive to the overall multiplicative bias in the shear 
measurements, we adopt the “conservative additive” selec¬ 
tion; this results in small additive systematic uncertainties, but 
possibly some moderate multiplicative systematic uncertain¬ 
ties. Eor ngmix, this selection removes galaxies with S/N<20 
and very small galaxies (Gaussian sigma smaller than the 
pixel scale). Eor imSshape, it removes galaxies with S/N< 15. 
In both cases, there were many other selections applied to both 
catalogs to remove stars, spurious detections, poor measure¬ 
ments, and various other effects that significantly biased shear 
estimates for both catalogs. 


3. The red-sequence Matched filter Galaxy Catalog (Redmagic) 

We use the DES SV red-sequence Matched-filter Galaxy 
Catalog (Redmagic Rozo et al., in preparation) v6.3.3 in this 
paper as one of the foreground samples. The objects in this 
catalog are photometrically selected luminous red galaxies 
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(LRGs). We use the terms Redmagic galaxies and LRG in¬ 
terchangeably. Specifically, Redmagic uses the Redmapper- 
calibrated model for the color of red-sequence galaxies as a 
function of magnitude and redshift [71]. This model is used 
to find the best-fit photometric redshift for all galaxies irre¬ 
spective of type, and the goodness-of-fit of the model is 
computed. For each redshift slice, all galaxies fainter than 
some minimum luminosity threshold Lmin are rejected. In ad¬ 
dition, Redmagic applies a selection where the 

^ function of redshift is chosen to ensure that the re¬ 
sulting galaxy sample has a nearly constant space density n. 
In this work, we set n = 10“^/z^Mpc“^. We assume fiat A 
CDM model with cosmological parameters = 0.7, /z = 100 
(varying these parameters does not change the results signif¬ 
icantly). The luminosity selection is L > 0.5L^(z), where the 
value of L^(z) at z=0.1 is set to match the Redmapper defi¬ 
nition for SDSS, and the redshift evolution for L^(z) is that 
predicted using a simple passive evolution starburst model 
at z = 3 [60]. We use the Redmagic sample because of the 
exquisite photometric redshifts of the Redmagic galaxy cat¬ 
alog: Redmagic photometric redshifts are nearly unbiased, 
with a scatter (Jz/(1 +z) ~ 1.7%, and a « 1.7% 4(7 redshift 
outlier rate. We refer the reader to Rozo et al. (in preparation) 
for further details of this catalog. 


B. Foreground and background galaxy samples selection 

As described in Sec. I, the main goal of this paper is to 
construct a projected mass map at a given redshift via weak 
lensing and to show that the mass map corresponds to real 
structures, or mass, in the foreground line-of-sight. For that 
purpose, we define two galaxy samples in this study — the 
background “source” sample which is lensed by foreground 
mass, and the foreground “lens” sample that traces the fore¬ 
ground mass responsible for the lensing. We wish to con¬ 
struct a weak lensing mass map from the background sample 
according to Sec. IIB and compare it with the mass map gen¬ 
erated from the foreground galaxy density map according to 
Sec. lie. 

We choose to have the two samples separated at redshift 
^ 0.55 in order to have a sufficient number of galaxies in both 
samples. Given that the photo-z training sample of our photo- 
z catalog does not extend to the same redshift and magnitude 
range as our data, we exclude objects with photo-z outside 
the range 0.1-1.2. The final foreground and background sam¬ 
ple are separated by the photo-z selection of 0.1 < z < 0.5 
and 0.6 < z < 1.2. Note that the Redmagic foreground galaxy 
sample has an additional redshift threshold z > 0.2. 

The main quantity of interest for the background galaxy 
sample is the shear measured for each galaxy. Since the back¬ 
ground sample only serves as a “backlight” for the foreground 
structure we are interested in, it need not be complete. There¬ 
fore the most important selection criteria for the background 
sample is to use galaxies with accurate shear measurements. 
Our source selection criteria are based on extensive tests of 
shear catalog as described in Jarvis et al. (in preparation). Af¬ 
ter applying the conservative selection of background galax¬ 


ies and our background redshift selection we are left with 
1,111,487 galaxies (2.22/arcmin^) for ngmix and 1,013,317 
galaxies (2.03/arcmin^) for imSshape. 




FIG. 1: Distributions of the single-point photo-z estimates for the 
background and foreground samples used in this paper are shown 
in the top panel, overlaid by the lensing efficiency (Eq. (10)) corre¬ 
sponding to the background sample. The background and the fore¬ 
ground main sample uses the mean of the PDF from BPZ for single¬ 
point estimates, while the LRG redshift estimate comes indepen¬ 
dently from Redmagic (see Sec. Ill A 3). The bottom panel shows the 
corresponding n(z) of the background and foreground main sample 
given by BPZ. These come from the sum of the PDF for all galaxies 
in the samples. 

The foreground sample in this work serves as the tracer of 
mass. Thus it is important to construct a magnitude-limited 
sample for which the number density is affected as little as 
possible by external factors. The main physical factors that 
contribute to variation in the galaxy number density are the 
spatial variation in depth and seeing. Both effects can intro¬ 
duce spatial variation in the foreground galaxy number den¬ 
sity, which can be wrongly identified as foreground mass fluc¬ 
tuations. We test both effects in Appendix A. Two subsam¬ 
ples are used in this work as foreground samples: the “main” 
foreground sample and the LRG foreground sample. While 
the space density of LRGs is significantly lower than that of 
the main sample, they are better tracers of galaxy clusters and 
groups, so we use them to check our results. The main fore¬ 
ground sample includes all the galaxies with i < 22 and the 
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TABLE I: Catalogs and selection criteron used to construct the foreground and background sample for this work, and the number of galaxies 
in each sample after all the selections are applied. All catalogs are based on the DES SVAl dataset. We use the Source Extractor MAG_AUTO 
parameter for the /-band magnitude. 



Background 

Foreground main 

Foreground LRG 

Input catalog 

ngmixOll imSshape 

SVAl Gold 

Redmagic 

Photometric redshift 

0.6<z<1.2 

0.1<z<0.5 

0.2<z<0.5 

Others 

“conservative additive” 

i<22 

constant density 
10-3 {h-^Mpc)-^ 

Number of galaxies 

1,111,487 1,013,317 

1,106,189 

28,033 

Number density (arcmin“^) 

2.22 

2.03 

2.21 

0.056 

Mean redshift 

0.826 

0.825 

0.367 

0.385 


LRG sample includes the LRGs in the Redmagic LRG cata¬ 
log with i <22. This magnitude selection is based on tests de¬ 
scribed in Appendix A 1 to ensure that our sample is shallower 
than the limiting magnitude for all regions of sky under study. 
The final main foreground sample contains 1,106,189 galaxies 
(2.21/arcmin^), while the LRG sample contains 28,033 galax¬ 
ies (0.05/arcmin^). Table I summarizes all the selection crite¬ 
ria applied on the three main samples used in this work. 

Fig. (1) shows the distributions of the single-point photo- 
z estimates (Zmean) for the final foreground and background 
samples overlaid by the lensing efficiency corresponding to 
the background sample (top panel), and the n{z) (from the 
BPZ code) for the background and main foreground sample 
(bottom panel). Note that the background galaxy number den¬ 
sity is much lower than the number density of all galaxies 
in the ngmixOll and imSshape catalogs, as we have made 
stringent redshift selections to avoid overlap between the fore¬ 
ground and background samples. 


C. Mock catalogs from simulations 

We use the simulations primarily as a tool to understand the 
impact of various effects on the expected signal, and a sanity 
check to confirm that our measurement method is producing 
reasonable results. We use a set of simulated galaxy catalogs 
“Aardvark vl.Oc” developed for the DES collaboration [72]. 
The full catalog covers approximately 1/4 of the sky and is 
complete to the final expected DES depth. 

The heart of the galaxy catalog generation is the algorithm 
Adding Density Determined Galaxies to Lightcone Simula¬ 
tions [ADDGALS; 72], which aims at generating a galaxy 
catalog that matches the luminosities, colors, and clustering 
properties of the observed data. The simulated galaxy catalog 
is based on three fiat ACDM dark matter-only N-body simula¬ 
tions, one each of a 1050 Mpc/h, 2600 Mpc/h and 4000 Mpc/h 
boxes with 1400^, 2048^ and 2048^ particles respectively. 
These boxes were run with LGadget-2 [73] with 2LPTic ini¬ 
tial conditions from [74] and GAME [75]. From an input lu¬ 
minosity function, galaxies are drawn and then assigned to 
a position in the dark matter simulation volume according to 
a statistical prescription of the relation between the galaxy’s 
magnitude, redshift and local dark matter density. The pre¬ 


scription is derived from a high-resolution simulation using 
SubHalo Abundance Matching techniques [72, 76, 77]. Next, 
photometric properties are assigned to each galaxy, where 
the magnitude-color-redshift distribution is designed to repro¬ 
duce the observed distribution of SDSS DR 8 [78] and DEEP2 
[79] data. The size distribution of the galaxies is magnitude- 
dependent and modelled from a set of deep (/ ^26) Suprime- 
Cam /-band images, which were taken at with seeing condi¬ 
tions of 0.6” [80]. Einally, the weak lensing parameters (k 
and 7 ) in the simulations are based on the ray-tracing al¬ 
gorithm Curved-sky grAvitational Lensing for Cosmological 
Light conE simulatioNS [CALCLENS; 81]. The ray-tracing 
resolution is accurate to 6.4 arcseconds, sufficient for the 
usage in this work. 

Aside from the intrinsic uncertainties in the modelling in 
the mock galaxy catalog (related to the input parameters and 
uncertainty in the galaxy-halo connection), there are also 
many real-world effects that are not included in these simula¬ 
tions, including as depth variation, seeing variation and shear 
measurement uncertainties. 


IV. MASS MAPS 

In Eig. (2) we show our final convergence maps generated 
using the data described in Sec. Ill A and the methods de¬ 
scribed in Sec. IIB and Sec. IIC. Eor the purpose of visu¬ 
alization we present maps for 20 arcmin Gaussian smoothing. 
In the top left panel we show the E-mode convergence map 
generated from shear. The top right panel shows the weighted 
foreground galaxy map from the main sample, Kg^main map. 
In both of these panels, red areas correspond to overdensities 
and blue areas correspond to under densities. The bottom left 
and bottom right panels show the product of the Ke (left) and 
Kb (right) maps with the Kg^main- Visually we see that there 
are more positive (correlated) areas for the Ke map compared 
to the Kb map, indicating clear detection of the weak lensing 
signal in these maps. Note that these positive regions could be 
either mass over-densities or under-densities. In Sec. VI, we 
present a quantitative analysis of this correlation. 

To estimate the significance of the structures in the mass 
maps, it is important to understand the noise properties of 
these maps. Uncertainties in the lensing convergence map 
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FIG. 2: The upper left panel shows the E-modes of the weak tensing convergence map. The upper right shows the weighted foreground galaxy 
map from the main sample, or Kg^^iain- The lower two panels show the product maps of the E-mode (left) and B-mode (right) convergence 
map with the Kg^rnain map. All maps are generated with a 5 arcmin pixel scale and 20 arcmin Gaussian smoothing. Red areas corresponds to 
overdensities and blue areas to underdensities in the upper panels. White regions correspond to the survey mask. The scale of the Gaussian 
smoothing kernel is indicated by the Gaussian profile on the upper right corner. 
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FIG. 3: The top panel shows the S/N map for the mass map in Fig. (2) estimated via randomized errors. Note that due to the Gaussian 
smoothing kernel, there is some mixing of scales which leads to higher contrasts in the cores of over and under-dense regions compared to 
top-hat smoothing. The bottom panel shows the normalized S/N distributions for both maps, overlaid by those measured from simulations 
described in Sec. VIB. The red dashed lines in both bottom panels show a Gaussian fit to the B-mode S/N. 


include contributions from both shape noise and measure¬ 
ment uncertainties, which is affected by the number density of 
galaxies across the field and the shear measurement method. 

We estimate the uncertainties on each pixel by randomising 
the shear measurements on each galaxy. A thousand random 
background galaxy catalogs were generated by shuffling the 
shear values between all the galaxies. We then construct Ke 
and Kb maps from these randomized catalogs in the same way 
as in Fig. (2). The standard deviation map for these 1000 ran¬ 
dom samples is used as the noise map. Dividing the signal 
map (Fig. (2)) by the noise map gives an estimate for the S/N 
of the different structures in the maps, as shown in Fig. (3). 


These values are broadly consistent with those predicted via 
Eq. (13) and simulations described in Sec. VIB. The bottom 
panels of Fig. (3) show the distribution of the S/N values for 
both E and B-mode maps for data as well as simulations pre¬ 
dicted by Eq. (13). We And that the B-mode distribution is 
consistent with a Gaussian distribution of standard deviation 
~ 1 as expected [82], and the E-mode gives more extreme 
values. The difference between the data and the simulation is 
consistent with cosmic variance and shape noise. 
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FIG. 4: The DBS SV mass map along with foreground galaxy clusters detected using the Redmapper algorithm. The clusters are overlaid as 
black circles with the size of the circles indicating the richness of the cluster. Only clusters with richness greater than 20 and redshift between 
0.1 and 0.5 are shown in the figure. The upper right corner shows the correspondence of the optical richness to the size of the circle in the plot. 
It can be seen that there is significant correlation between the mass map and the distribution of galaxy clusters. Several superclusters (black 
squares) and voids (white squares) can be identified in the joint map. 


V. CORRELATION WITH GALAXY CLUSTERS AND 
POTENTIAL SUPER STRUCTURES 

In this section we compare our mass map with optically 
identified groups and clusters of galaxies using the Redmap¬ 
per algorithm (Rykoff et al. in preparation) on DBS data. We 
overlay in Fig. (4) Redmapper clusters and groups on the mass 
map as black circles. The size of these circles corresponds to 
the optical richness of these structures. Only clusters with 
optical richness A greater than 20 and redshift between 0.1 
and 0.5 are shown in the figure. According to Rykoff et al. 
[83] and Saro et al. (in preparation), cluster mass scales 


approximately linear with A, with A = 20 corresponding to 
^ 1.7 X 10^"^ Mq and A = 80 corresponding to ^ 7.6 x 10^"^ 
Mq. It is evident from this figure that the structures in the 
weak lensing mass map have significant correlation with the 
distribution of optically identified Redmapper clusters. The 
combination of the lensing mass maps, Redmapper clusters 
and Redmagic LRGs provides a powerful tool for identifying 
superstructures in the universe that would otherwise be hard 
to spot. 

Superclusters are the largest distinct structures in the uni¬ 
verse, typically 10 Mpc or larger in extent with fractional 
overdensities of order 1-10 times the cosmic mean density. 
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FIG. 5: Left: the blue curve in each of the 4 panels shows the weighted redshift distribution of galaxy clusters counts with optical richness 
A >5 at the 4 different locations in the mass map of Fig. (4) corresponding to large convergence peak locations. The RA, Dec coordinates of 
these pointings are shown in the top right corner of each panel and the field numbers are listed on the top left corner. The counts are calculated 
for a 1 deg radius area, and the histograms are weighted by A and the lensing efficiency to properly represent the mass distribution and the 
lensing probed by the mass map. The thick grey line indicates the corresponding average number count in the full map. The redshift range 
above z = 0.6 is marked with the shaded grey area, as these ranges overlap with the background sample. Right: a candidate supercluster is 
shown by zooming in on a narrow redshift range of field 2 (red band in upper right panel on the left) where a peak in the cluster counts occurs. 
Each circle indicates the location of a galaxy belonging to a Redmapper cluster. The large spatial extent (a transverse distance of 10 Mpc is 
indicated in the panel) and the irregular shape characteristic of 3D superclusters is evident. 
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FIG. 6: Left: same as the left panel of Fig. (5) but plotted for voids identified in the mass map. There are typically fewer than average clusters 
over much of the line of sight which also contains some deep underdense regions at specific redshifts. At the higher redshifts, there are also 
above average cluster counts, but since the redshift range overlaps with the source galaxy sample, the interpretation of the structures is more 
complicated. Right: radial distribution of the Redmagic LRGs for field 5 in the left panel (red bands in upper left panel). The data are consistent 
with the existence of two voids modeled by the “top-hat” void model of width 190 Mpc//z and 120 Mpc//z respectively. 


Cosmic voids are the corresponding underdensities, typically 
larger than 10 Mpc in radius with fractional underdensity of 
order unity. We identify superclusters and voids from the mass 
and galaxy maps in Fig. (2) and Fig. (4). The large peaks at 
the positions (RA, Dec) = (71.0, -45.0), (69.9, -47.8), (69.7, 
-54.5) and (69.1, -57.3) and large voids at (RA, Dec) = (65.6, - 


49.0), (75.1, -54.6), (75.7, -58.0) and (82.8, -59.5) are selected 
as shown in Fig. (4). The transverse spatial extent of these su¬ 
perstructures is typically greater than 10 Mpc. We compare 
in the left panels of Fig. (5) and Fig. (6) the redshift distribu¬ 
tion of the foreground clusters within 1 deg radius of these 
locations with the average redshift distribution of the clus- 
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ters in the entire SV field. The histograms are weighted by 
the optical richness A as well as the tensing efficiency of our 
source sample (Fig. (1)). A scales roughly linearly with the 
total mass of the cluster [83]. We find that some of the mass 
map peaks correspond to supercluster-like structures that are 
localized in narrow redshift ranges, while others (e.g. field 
3) show evidence for structures extending over wider redshift 
range. On the other hand the large voids typically have fewer 
clusters than average along the line of sight and some deep 
underdense regions (candidate 3D voids) at specific redshifts. 
In some cases there are also above average cluster counts in 
small ranges in redshift (field 6), as expected from the pro¬ 
jected nature of these mass maps. The redshift range above 
z = 0.6 is marked with the shaded grey area, as this range 
overlaps with the background sample thus complicating the 
interpretation of the relationship with the mass map. In the fu¬ 
ture we will carry out more detailed studies of the mass maps 
using lensing tomography. 

We show two cases for further investigations of potential 
superclusters and voids identified through this method. First, 
we look at the spatial distribution of the cluster members in 
thin redshift slices, identical to the analysis in Melchior et al. 
[70], and find structures such as the one shown in the right 
panel of Fig. (5). The redshift extent Az =0.03 corresponds to 
a line-of-sight distance of about 90 Mpc//z, while the trans¬ 
verse size of the structure shown on the right is about 20 
Mpc//z. The line of sight scale corresponds to the size of 
the largest filamentary structures in cosmological simulations 
[84]. These numbers indicate that this is a good candidate for 
a 3D supercluster. The tight photo-z accuracy of the Redmap- 
per clusters ^ 0.01 (1 +z)) gives us confidence in the iden¬ 
tification of real 3D structures. 

For the voids, we follow the method developed in Szapudi 
et al. [85] and study the radial distribution of the foreground 
Redmagic LRGs. We use LRGs within 0.5 deg radius of the 
chosen position and calculate 5lrg = {^lrg — ^lrc^I^lrg in 
100 Mpc//z radial bins, where ulrg is the number of LRGs in 
that bin and hlrg is the average number of LRGs for the full 
Redmagic catalog in the same radial bin. The radial profile for 
one void is shown in the right panel of Fig. (6): it is consistent 
with two large voids in this line of sight. We use a simple 
“top-hat” void model [85] with an amplitude 5lrg = —0.7, 
an extent of 190 Mpc/h at a distance of 750 Mpc/h for the 
first void, and another one with Srrg = —0.7, an extent of 120 
Mpc/h at 1250 Mpc/h. The combination of these two void 
models, smoothed by the photo-z uncertainty, matches well 
with the data. We also observe that there could be a similarly 
large but shallower void at higher redshift, also contributing 
to the projected underdensity in the mass map. 

The size and mass of the superclusters are of interest for 
cosmology as they represent the most massive end of the 
matter distribution. The is especially interesting as the DBS 
dataset allows us to extend our studies to z ~ L We defer 
more detailed studies of superclusters and voids to follow up 
work. 


VI. CORRELATION WITH GALAXY DISTRIBUTION 

In this section we quantitatively analyze the extent to which 
mass follows galaxy density in the data. To do this, we cross¬ 
correlate the weak lensing mass map with the weighted fore¬ 
ground galaxy density map. The correlation is quantified 
via the Pearson cross-correlation coefficient as described in 
Sec. VIA. We cross check the results using simulations in 
Sec. VIB. 


A. Quantifying the galaxy-mass correlation 


We smooth both the convergence maps generated from 
weak lensing and from the foreground galaxy density with a 
Gaussian filter. These smoothed maps are used to estimate 
the correlation between the foreground structure and the weak 
lensing convergence maps. We calculate the correlation as a 
function of the smoothing scale. The correlation is quantified 
via the Pearson correlation coefficient defined as 


pKEKg — 


i^EKg) 
(^Ke ’ 


( 16 ) 


where {^RKg) is the covariance between Kr and Kg', g^e and 
G^^ are the standard deviation of the Kr map, and the Kg map 
from either the foreground main galaxy sample or the fore¬ 
ground LRG sample. In this calculation, pixels in the masked 
region are not used. We also remove pixels within 10 ar- 
cmin of the boundaries to avoid significant artefacts from the 
smoothing. 

Fig. (7) shows the Pearson correlation coefficient as func¬ 
tion of smoothing scales from 5 to 40 arcmin. We find that 
there is significant correlation between the weak lensing E- 
mode convergence and convergence from different foreground 
samples, with increasing correlation towards large smooth¬ 
ing scale. This trend is expected for noise-dominated maps, 
because the larger smoothing scales reduce the noise fiuc- 
tuations in the map significantly. A similar trend is found 
when using the LRGs as foreground instead of the general 
magnitude-limited galaxy sample. The lower Pearson corre¬ 
lation between the mass map and LRG sample is because of 
the larger shot noise due to the lower number density com¬ 
pared to the magnitude-limited foreground sample. The error 
bar on the correlation coefficient is estimated based on jack¬ 
knife resampling. We divide the observed sky into jackknife 
regions of size 10 deg^ and recalculate the Pearson correla¬ 
tion coefficients, excluding one of the 10 deg^ regions each 
time. We found that the estimated uncertainties do not de¬ 
pend significantly on the exact value of patch size. We esti¬ 
mate the correlation coefficient after removing one of those 
patches from the sample to get jackknife realizations of the 
cross-correlation coefficient pj. Finally, the variance is esti¬ 
mated as 


^P = ^-R^'L^Pj-P?^ ( 17 ) 


where j runs over all the V jackknife realizations and p is the 
average correlation coefficients of all patches. 
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FIG. 7: This figure shows the Pearson correlation coefficient between foreground galaxies and convergence maps as a function of smoothing 
scale for the ngmix galaxy catalog. The solid and open symbols show the E and B-mode correlation coefficients respectively. The black 
circles are for the main foreground sample and the red circles for foreground LRGs. The grey shaded regions show the la bounds for E and 
B mode correlations from simulations for the main foreground sample with the same pixelization and smoothing (see Sec. VIB for details). 
We do not show the similar simulation results for the LRG sample. The detection significance for the correlation is in the range ^ 5 — 7(7 at 
different smoothing scales. The green points show the correlation between E and B-modes of the mass map. The various B-mode correlations 
are consistent with zero. Uncertainties on all measurements are estimated based on jackknife resampling. 



EIG. 8: Same as Eig. (7) but using the imSshape galaxy catalog. 


We find that the Pearson correlation coefficient between 
Kg from the main foreground galaxy sample (LRG sample) 
and weak lensing E-mode convergence is 0.39 ± 0.06 (0.36 ± 
0.05) at 10 arcmin smoothing and 0.52 ±0.08 (0.46 ±0.07) at 
20 arcmin smoothing. This corresponds toa^6.8(7 (7.5 a) 
significance at 10 arcmin smoothing and ~6.8(7(6.4(7)at20 
arcmin smoothing. As a zeroth-order test of systematics we 
also estimated the correlation between the B-mode weak lens¬ 


ing convergence and the Kg maps. We find that the correla¬ 
tion between k^ and the main foreground sample is consistent 
with zero at all smoothing scales. Similarly, the correlation 
between E and B modes of k is consistent with zero. Eor com¬ 
parison, we show the same plot calculated for the imSshape 
catalog in Eig. (8). We find very similar results, with slightly 
larger correlation between Ke and Kb at the I a level. 
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FIG. 9: Maps from simulations that are designed to mimic the data 
in our analysis. The simulations are generated for a field of size 
15 X 17.6 deg^ with similar redshift and magnitude selections for the 
foreground and the background sample as the data. The true K and 
Kg maps are shown in the first row, where Kg is modelled for the main 
foreground sample. The reconstructed Ke and Kb maps from the true 
7 are shown in the first two panels of the second row, followed by the 
Ke and Kb maps reconstructed from the ellipticity (e) values. The 
last row first shows the Ke and Kb constructed from e with photo-z 
uncertainties, then the same maps with an SV survey mask applied. 
The last two panels on the bottom most closely match the data. 
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FIG. 10: Pearson correlation coefficient pxKg between the different 
simulated maps shown in Fig. (9) as a function of smoothing scale. X 
represents the different K maps as listed in the legend. This plot is the 
simulation version of Fig. (7), where one can see how the measured 
values in the data could have been degraded due to various effects. 
The qualitative trend of the correlation coefficients as a function of 
smoothing scale is consistent with that observed in data. When re¬ 
constructing Ke from the true 7 small errors are introduced due to the 
nonlocal reconstruction, lowering the correlation coefficient by a few 
percent. Adding shape noise to the shear measurement lowers the 
signal significantly, with the level of degradation dependent on the 
smoothing scale. Adding photo-z uncertainties changes the signal by 
a few percent. Finally, placing an SV-like survey mask changes the 
signal by ~10%. The black curve with its error bars corresponds to 
the shaded region in Fig. (7). 


B. Comparison with mock catalogs 

At this point, it is important to verify whether our mea¬ 
surements in the data are consistent with what is expected. 
We investigate this using the simulated catalogs described in 
Sec. Ill C. As the simulations lack several realistic systematic 
effects in the data, these tests mainly serve as a guidance for 
us to understand: (1) the origin of the B-mode in the K maps, 
(2) the approximate expected level of p^Xg under pixelization 
and smoothing, (3) the effect on p^Kg from photo-z uncertain¬ 
ties and cosmic variance, and (4) the effect on the maps and 
pKKg from the survey mask. 

We construct a sample similar to the SV data. The same 
redshift, magnitude, and number density selections in Table I 
are applied to the simulations to form a foreground and a back¬ 
ground sample. We choose to simulate the main foreground 
sample as the LRG foreground sample selection in the sim¬ 
ulations is less controllable. For the background sample, we 
add Gaussian noise with standard deviation a =0.27 to each 
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component of the true shear in the simulations to generate a 
model for the ellipticities that matches the data (Jarvis et al. 
in preparation). We then create a Kg map from the main fore¬ 
ground sample and a k map from the background sample the 
same way as is done in the data. The cross-correlation co¬ 
efficient pxKg is calculated from these simulated maps as in 
Sec. VI. We consider the same range of smoothing scales for 
the maps when calculating as that in Fig. ( 8 ). 

The simulations provide us a controlled way of separating 
the different sources of effects. We construct the maps in the 
following steps, in order of increasing similarities to data: ( 1 ) 
pixelating and smoothing the true K values; (2) constructing 
the K values from the true 7 values; (3) construct the k values 
from the galaxy ellipticities which include shape noise (we 
generate 20 realizations); (4) using a photo-z model for the 
foreground and the background galaxies instead of the true 
redshift; generate four different maps from different regions 
on the sky; ( 6 ) use the SV survey mask. Note that in step (3) 
we take the galaxy ellipticity to be the sum of a random com¬ 
ponent (sampled from a Gaussian with standard deviation of 
0.27) and the lensing shear, this model is designed to match 
the data, which includes the intrinsic shape noise and other 
measurement noise associated with e.g. the PSF modelling. In 
step (4) we have modelled the photo-z errors from a spectro¬ 
scopic sample that ran through the same photo-z code, taking 
the spectroscopic redshift to be the “true” redshift. 

The difference between step (1) and step (2) measures the 
quality of the KS reconstruction method. The difference be¬ 
tween step (2) and step (3) shows the effect of shape noise and 
measurement noise. Steps (4), (5) and ( 6 ) then show the effect 
of photo-z uncertainties, cosmic variance and masking. For 
each SV-size maps, we generate 20 (shape noise) x 4 (cosmic 
variance) X 2 (photo-z) x2 (mask)=320 corresponding simu¬ 
lations. 


1. Maps from simulations 

Fig. (9) shows the various maps generated from one partic¬ 
ular patch of the simulations in this procedure for 5 arcmin 
pixels and 20 arcmin smoothing scales (consistent with that 
in Fig. (2)). The amplitude of Kp and Kp both become larger 
than in the true maps when shape noise is added, and the re¬ 
sulting Kp map has only slightly higher contrast than the Kp 
map. When photo-z uncertainties are included, we see that the 
peaks and voids in the Kp maps visibly move around. Apply¬ 
ing the mask mainly changes the morphology of the structures 
in the maps around the edges. Comparing the last Kp panel in 
Fig. (9) and Fig. (2), we see that the amplitude and qualitative 
scales of the variation in the Kp maps are similar. On the other 
hand, if we compare the Kg maps in the simulations with the 
Kg maps in Fig. (2), we find some qualitative differences be¬ 
tween the simulations and the data. The simulation contains 
more small scale structure and low-x*^ regions compared to 
the data. We do not investigate this issue further here, as the 
level of agreement in the simulations and the data is sufficient 
for our purpose. 


2. Correlation coefficients from simulations 

Fig. (10) shows the mean Pearson correlation coefficient 
between the different maps as a function of smoothing scales 
for the 80 sets of simulated maps (4 different areas in the sky 
and 20 realisations of shape noise each). The error bars indi¬ 
cate the standard deviation of these 80 simulations. 

We find pKtrue^^g ~ 0-8 — 0.9. Several factors contribute to 
this. First, the foreground galaxy sample only includes a finite 
redshift range, and not all galaxies that contribute to the Ktrue 
map. Second, the presence of a redshift-dependent galaxy bias 
adds further complication to the correlation coefficient. The 
effect of converting from the true shear 7 to convergence low¬ 
ers the correlation coefficient by about 3%. This is a measure 
of the error in the KS conversion under finite area and reso¬ 
lution of the shear fields. The main degradation of the signal 
comes when shape noise and measurement noise is included. 
Photo-z uncertainties in both the foreground and the back¬ 
ground sample changes the correlation coefficient slightly. Fi¬ 
nally, the survey mask lowers the correlation coefficient by 
- 10 %. 

The final correlation coefficient after considering all the ef¬ 
fects discussed above is shown by the black curve in Fig. (10) 
and overplotted as the shaded region in Fig. (7). We find that 
the dependence of p^Kg on the smoothing scale in the simula¬ 
tion is qualitatively and quantitatively very similar to that seen 
in Fig. (7). 


VII. SYSTEMATIC EFFECTS 

In this section we examine the possible systematic uncer¬ 
tainties in our measurement. We focus on the cross correla¬ 
tion between our weak lensing mass map Kp and the main 
foreground density map Kg^main- To simplify the notation, 
we omit the “main” in the subscript and use Kg to represent 
the main foreground map in this section. We investigate the 
potential contamination from systematic effects on the cross¬ 
correlation coefficient pKEKg by looking at the spatial correla¬ 
tion of various quantities with the Kp map and the Kg map. 

As discussed in Appendix A, there are several factors that 
can contaminate the 5g maps. For example, depth and PSF 
variations in the observed field can introduce artificial clus¬ 
tering in the foreground galaxy density map. Although we 
use magnitude and redshift selections according to the tests in 
Appendix A, one can expect some level of residual effects on 
the Kg maps. The Kp map is constructed from shear catalogs 
of the background sample, thus systematics in the shear mea¬ 
surement will propagate into the Kp map. In Jarvis et al. (in 
preparation), extensive tests of systematics have been carried 
out on the shear catalog. Therefore here we focus on the sys¬ 
tematics that are specifically relevant for mass mapping and 
the correlation coefficient Eq. (16). 

We identify several possible sources of systematics for 
the background and foreground sample as listed in Table II. 
We generate maps of these quantities that are pixelated and 
smoothed on the same scale as the Kp and Kg maps. We then 
evaluate the contribution of these effects to the correlation co- 
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FIG. 11: The normalised cross-correlation coefficient is shown for 20 different systematic uncertainty parameters. The systematics 

parameters, represented by 0, are listed in Table II and shown for two smoothing scales. The Pk-^k-^;© values are normalized by Pk-^k-^ to show 
the relative magnitude of the systematic and the signal. The red dashed line indicates where the systematic is 5% of PxEXg- The error bars are 
estimated from resampling the foreground and background galaxy sample in patches of size 10 deg^. The left panel is calculated for ngmix 
while the right panel is for im3shape. 



FIG. 12: Pearson correlation coefficient Pk-^© where 0 represents the quantities listed in Table II. We show the statistics for two smoothing 
scales and for both ngmix (left) and imSshape (right). The right-most points in both panel correspond to the detection signal in Fig. (7) and 
Fig. (8). The error bars are estimated from resampling the foreground and background galaxy sample in patches of size 10 deg^. Note that this 
is a different statistic from that in Fig. (11), thus the y-axis values are not directly comparable. 


efficient (Eq. (16)) based on the following diagnostic quantity: 


P/C£/c„;0 — 


PKE&PKg& 

pee 


(18) 


with pxY being the cross-correlation function, which is ef¬ 
fectively the unnormalized Pearson correlation coefficient be¬ 
tween X and Y, or 


Pxy = {XY). (19) 


Equation 18 measures the contribution from some system¬ 


atics field 0 to pKEKg- We calculate pKEKg-,e with 0 being 
any of the 20 quantities in Table II (excluding the signal). 
Eig. (11) shows the normalized cross-correlation coefficient 
pKEKg\e/pKEKg values for all the quantities considered for 10 
and 20 arcmin smoothing, with the red dashed line at 5%. The 
error bars are estimated by jackknife resampling similar to that 
described in Sec. VIA, and the two panels show the results 
for ngmix and imSshape respectively. The normalized cross¬ 
correlation coefficient is a measure of the fractional contam¬ 
ination in the Pearson coefficient (Eq. (16)) from each of the 
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TABLE II: Quantities examined in our systematics tests. 


Map name 

Description 

kE (signal) Ke from ji , 72 for background sample 

kg (signal) 

Kg from main foreground sample 

kB 

Kb from 71 , 72 for background sample 

ns_f 

star number per pixel 

ng-b 

galaxy number per pixel for background sample 

snr 

signal-to-noise of galaxies in imSshape 

mask 

fraction of area masked in galaxy postage stamp 

gl 

average 71 for background sample 

g 2 

average 72 for background sample 

psf_el 

average PSF ellipticity 

psf_e 2 

average PSF ellipticity 

psf_T 

average PSF size (ngmix only) 

psfTwhm 

average PSF size (imSshape only) 

psf_kE 

Ke generated from average PSF ellipticity 

psfJcB 

Kb generated from average PSF ellipticity 

zp_b 

mean photo-z for background sample 

zp± 

mean photo-z for foreground sample 

ebv 

mean extinction 

sky sigma 

standard deviation of sky brightness in ADU 

sky 

mean sky brightness in ADU 

maglim 

mean limiting /-band AB magnitude 

exptime 

mean exposure time in seconds 

airmass 

mean airmass 


systematics maps 0. 

We find that for ngmix all quantities show contributions to 
the systematic uncertainties at 10 arcmin smoothing to be at 
the level of 5% or lower, while the systematics increase to 
up to 15% when smoothing at the 20 arcmin scale (though 
with large error bars on the systematics estimation). For 
imSshape, most of the values stay below 5% for both smooth¬ 
ing scales. The largest contribution in both cases come from 
the variation in the PSF properties (psLel, psf_e2, psLkB). 
This is expected, as the modelling of the PSF is known to be 
a significant challenge in weak lensing. Since all these PSF 
quantities are correlated with each other, and many other pa¬ 
rameters (gl, g2, snr, maglim) are correlated with the PSF 
properties, we do not expect the total systematics contami¬ 
nation to be a direct sum of all these parameters. Instead, 
we discuss in Appendix B how one can isolate the indepen¬ 
dent contributions of the systematics via a Principal Compo¬ 
nent Analysis approach and correct for them. We find that the 
correction changes the final Pearson correlation coefficient by 
3.5% relative to the original pxEXg measured in Sec. VI. 

Finally, to check the level of systematic contamination in 
our Ke map itself, we also calculate the Pearson correlation 
coefficient (Eq. (16)) between the various maps in Table II and 
our Ke map. Note that this contamination may or may not be 
pronounced in Fig. (11) since the statistics plotted there also 
take into account the correlation of Kg with the various quanti¬ 
ties. This test is independent of the foreground map, therefore 
is important for applications of the Ke map that do not also use 
the foreground maps. Fig. (12) shows the resulting 21 Pear¬ 
son correlation coefficients. We find that the signal shown in 
the right-most points in the plot (pKsKg) is larger than all other 


correlations by at least a factor of ~3. 

We also note that in both of these tests, the area of the map 
is not big enough to ignore the fact that some of these cor¬ 
relations can be intrinsically non-zero, even if there were no 
systematics contamination in the maps. 


VIII. CONCLUSIONS 

In this work, we present a weak lensing mass map based on 
galaxy shape measurements in the 139 deg^ SPT-E field from 
the Dark Energy Survey Science Verification data. We have 
cross-correlated the mass map with maps of galaxy and cluster 
samples in the same dataset. We demonstrate that candidate 
superclusters and voids along the line of sight can be identified 
exploiting the tight scatter of the cluster photo-z’s. 

We constructed mass maps from the foreground Redmagic 
LRG and general magnitude-limited galaxy samples under the 
assumption that mass traces light. We find that the E-mode of 
the convergence map correlates with the galaxy based maps 
with high statistical significance. We repeated this analysis 
at various smoothing scales and compared the results to mea¬ 
surements from mock catalogs that reproduce the galaxy dis¬ 
tribution and lensing shape noise properties of the data. The 
Pearson cross-correlation coefficient is 0.39 ± 0.06 (0.36 ± 
0.05) at 10 arcmin smoothing and 0.52 ±0.08 (0.46 ±0.07) 
at 20 arcmin smoothing for the main (LRG) foreground sam¬ 
ple. This corresponds toa^6.8(7 {1.5a) significance at 10 
arcmin smoothing and ^6.8(7 (6.4(7) at 20 arcmin smooth¬ 
ing. We get comparable values from the mock catalogs, indi¬ 
cating that statistical uncertainties, not systematics, dominate 
the noise in the data. The B-mode of the mass map is consis¬ 
tent with noise and its correlations with the foreground maps 
are consistent with zero at the 1(7 level. 

To examine potential systematic uncertainties in the conver¬ 
gence map we identified 20 possible systematic tracers such 
as seeing, depth, PSF ellipticity and photo-z uncertainties. We 
show that the systematics effects are consistent with zero at 
the 1 or 2(7 level. In Appendix B, we present a simple scheme 
for the estimation of systematic uncertainties using Principal 
Component Analysis. We discuss how these contributions can 
be subtracted from the mass maps if they are found to be sig¬ 
nificant. 

The results from this work open several new directions of 
study. Potential areas include the study of the relative distribu¬ 
tion of hot gas with respect to the total mass based on X-ray 
or SZ observations, estimation of galaxy bias, constraining 
cosmology using peak statistics, and finding filaments in the 
cosmic web. The tools that we have developed in this paper 
are useful both for identifying potential systematic errors and 
for cosmological applications. The observing seasons for the 
first two years of DBS are now complete [86] and survey an 
area well over ten times that of the SV data, though shallower 
by about half a magnitude. The full DBS survey area will 
be ^ 35 times larger than that presented here, at roughly the 
same depth. The techniques and tools developed in this work 
will be applied to this new survey data, allowing significant 
expansion of the work presented here. 



18 


Acknowledgements 

We are grateful for the extraordinary contributions of our 
CTIO colleagues and the DECam Construction, Commission¬ 
ing and Science Verification teams in achieving the excellent 
instrument and telescope conditions that have made this work 
possible. The success of this project also relies critically on 
the expertise and dedication of the DES Data Management 
group. 

We thank Jake VanderPlas, Andy Connolly, Phil Marshall, 
Ludo van Waerbeke, and Rafal Szepietowski for discussions 
and collaborative work on mass mapping methodology. CC 
and AA are supported by the Swiss National Science Eoun- 
dation grants 200021-149442 and 200021-143906. SB and 
JZ acknowledge support from a European Research Coun¬ 
cil Starting Grant with number 240672. DG was supported 
by SEB-Transregio 33 ‘The Dark Universe’ by the Deutsche 
Eorschungsgemeinschaft (DEG) and the DEG cluster of ex¬ 
cellence ‘Origin and Structure of the Universe’. ES acknowl¬ 
edges financial support provided by CAPES under contract 
No. 3171-13-2. OL acknowledges support from a European 
Research Council Advanced Grant EP7/291329 

Eunding for the DES Projects has been provided by the U.S. 
Department of Energy, the U.S. National Science Eoundation, 
the Ministry of Science and Education of Spain, the Science 
and Technology Eacilities Council of the United Kingdom, the 
Higher Education Eunding Council for England, the National 
Center for Supercomputing Applications at the University of 
Illinois at Urbana-Champaign, the Kavli Institute of Cosmo¬ 
logical Physics at the University of Chicago, the Center for 
Cosmology and Astro-Particle Physics at the Ohio State Uni¬ 
versity, the Mitchell Institute for Eundamental Physics and 
Astronomy at Texas A&M University, Einanciadora de Es- 


tudos e Projetos, Punda 9 ao Carlos Chagas Eilho de Amparo 
a Pesquisa do Estado do Rio de Janeiro, Conselho Nacional 
de Desenvolvimento Cientifico e Tecnologico and the Min- 
isterio da Ciencia e Tecnologia, the Deutsche Eorschungsge¬ 
meinschaft and the Collaborating Institutions in the Dark En¬ 
ergy Survey. 

The DES data management system is supported by the 
National Science Eoundation under Grant Number AST- 
1138766. The DES participants from Spanish institutions 
are partially supported by MINECO under grants AYA2012- 
39559, ESP2013-48274, EPA2013-47986, and Centro de Ex- 
celencia Severo Ochoa SEV-2012-0234, some of which in¬ 
clude ERDE funds from the European Union. 

The Collaborating Institutions are Argonne National Lab¬ 
oratory, the University of California at Santa Cruz, the Uni¬ 
versity of Cambridge, Centro de Investigaciones Energeticas, 
Medioambientales y Tecnologicas-Madrid, the University of 
Chicago, University College London, the DES-Brazil Con¬ 
sortium, the Eidgenossische Technische Hochschule (ETH) 
Zurich, Eermi National Accelerator Laboratory, the Uni¬ 
versity of Edinburgh, the University of Illinois at Urbana- 
Champaign, the Institut de Ciencies de I’Espai (lEEC/CSIC), 
the Institut de Eisica d’Altes Energies, Lawrence Berke¬ 
ley National Laboratory, the Ludwig-Maximilians Universitat 
and the associated Excellence Cluster Universe, the Univer¬ 
sity of Michigan, the National Optical Astronomy Observa¬ 
tory, the University of Nottingham, The Ohio State University, 
the University of Pennsylvania, the University of Portsmouth, 
SLAC National Accelerator Laboratory, Stanford University, 
the University of Sussex, and Texas A&M University. 

This paper has gone through internal review by the DES 
collaboration. 


[1] M. Bartelmann and P. Schneider, Physics Reports 340, 291 
(2001), astro-ph/9912508. 

[2] H. Hoekstra and B. Jain, Annual Review of Nuclear and Particle 
Science 58 , 99 (2008), 0805.0139. 

[3] L. Van Waerbeke, Y. Mellier, T. Erben, J. C. Cuillandre, 
K Bernardeau, R. Maoli, E. Bertin, H. J. McCracken, O. Le 
Fevre, B. Fort, et al., A&A 358 , 30 (2000), astro-ph/0002500. 

[4] N. Kaiser, G. Wilson, and G. A. Luppino, ArXiv Astrophysics 
e-prints (2000), astro-ph/0003338. 

[5] D. J. Bacon, A. R. Refregier, and R. S. Ellis, MNRAS 318, 625 
(2000), astro-ph/0003008. 

[6] H. Hoekstra, Y. Mellier, L. van Waerbeke, E. Semboloni, L. Fu, 
M. J. Hudson, L. C. Parker, I. Tereno, and K. Benabed, ApJ 
647 , 116 (2006), astro-ph/0511089. 

[7] H. Lin, S. Dodelson, H.-J. Seo, M. Soares-Santos, J. Annis, 
J. Hao, D. Johnston, J. M. Kubo, R. R. R. Reis, and M. Simet, 
ApJ 761 , 15 (2012), 1111.6622. 

[8] C. Heymans, L. Van Waerbeke, L. Miller, T. Erben, H. Hilde- 
brandt, H. Hoekstra, T. D. Kitching, Y. Mellier, P. Simon, 
C. Bonnett, et al., MNRAS 427, 146 (2012), 1210.0032. 

[9] D. J. Bacon, R. J. Massey, A. R. Refregier, and R. S. Ellis, 
MNRAS 344 , 673 (2003), astro-ph/0203134. 

[10] M. Jarvis, B. Jain, G. Bernstein, and D. Dolney, ApJ 644 , 71 


(2006), astro-ph/0502243. 

[11] E. Semboloni, Y. Mellier, L. van Waerbeke, H. Hoekstra, 
I. Tereno, K. Benabed, S. D. J. Gwyn, L. Fu, M. J. Hudson, 

R. Maoli, et al., A&A 452, 51 (2006), astro-ph/0511090. 

[12] L. Fu, M. Kilbinger, T. Erben, C. Heymans, H. Hildebrandt, 
H. Hoekstra, T. D. Kitching, Y. Mellier, L. Miller, E. Sem¬ 
boloni, et al., MNRAS 441 , 2725 (2014), 1404.5469. 

[13] M. J. Jee, J. A. Tyson, M. D. Schneider, D. Wittman, 

S. Schmidt, and S. Hilbert, ApJ 765, 74 (2013), 1210.2732. 

[14] M. Kilbinger, L. Fu, C. Heymans, F. Simpson, J. Benjamin, 

T. Erben, J. Harnois-Deraps, H. Hoekstra, H. Hildebrandt, T. D. 
Kitching, et al., MNRAS 430 , 2200 (2013), 1212.3338. 

[15] L. Van Waerbeke, G. Hinshaw, and N. Murray, Phys.Rev.D 89 , 
023508 (2014), 1310.5721. 

[16] B. Jain and L. Van Waerbeke, ApJ 530, LI (2000), astro- 
ph/9910459. 

[17] P. Fosalba, E. Gaztanaga, F. J. Castander, and M. Manera, MN¬ 
RAS 391 , 435 (2008), 0711.1540. 

[18] J. P. Dietrich and J. Hartlap, MNRAS 402, 1049 (2010), 
0906.3512. 

[19] J. M. Kratochvil, Z. Haiman, and M. May, Phys.Rev.D 81 , 
043519(2010), 0907.0486. 

[20] J. Berge, A. Amara, and A. Refregier, ApJ 712 , 992 (2010), 



19 


0909.0529. 

[21] A. Amara, S. Lilly, K. Kovac, J. Rhodes, R. Massey, 

G. Zamorani, C. M. Carollo, T. Contini, J.-R Kneib, O. Le 
Fevre, et al., MNRAS 424 , 553 (2012), 1205.1064. 

[22] M. Jauzac, E. Julio, J.-P. Kneib, H. Ebeling, A. Leauthaud, C.- 
J. Ma, M. Limousin, R. Massey, and J. Richard, MNRAS 426 , 
3369 (2012), 1208.4323. 

[23] H. Y. Shan, J.-P. Kneib, J. Comparat, E. Julio, A. Charbonnier, 
T. Erben, M. Makler, B. Moraes, L. Van Waerbeke, E. Courbin, 
et al., MNRAS 442 , 2534 (2014), 1311.1319. 

[24] H. S. Hwang, M. J. Geller, A. Diaferio, K. J. Rines, and H. J. 
Zahid, ApJ 797, 106 (2014), 1410.3883. 

[25] M. White, L. van Waerbeke, and J. Mackey, ApJ 575, 640 
(2002), astro-ph/0111490. 

[26] C. Heymans, M. E. Gray, C. Y. Peng, L. van Waerbeke, E. E. 
Bell, C. Wolf, D. Bacon, M. Balogh, E. D. Barazza, M. Barden, 
et al., MNRAS 385 , 1431 (2008), 0801.1156. 

[27] N. A. Bahcall and X. Ean, ApJ 504,1 (1998), astro-ph/9803277. 

[28] Z. Haiman, J. J. Mohr, and G. P. Holder, ApJ 553 , 545 (2001), 
astro-ph/0002336. 

[29] G. Holder, Z. Haiman, and J. J. Mohr, ApJ 560 , LI 11 (2001), 
astro-ph/0105396. 

[30] L. Knox, Y.-S. Song, and J. A. Tyson, Phys.Rev.D 74 , 023512 
(2006), astro-ph/0503644. 

[31] B. Jain and J. Khoury, Annals of Physics 325 , 1479 (2010), 
1004.3294. 

[32] M. Plionis and S. Basilakos, MNRAS 330, 399 (2002), astro- 
ph/0106491. 

[33] L. Van Waerbeke, J. Benjamin, T. Erben, C. Heymans, 

H. Hildebrandt, H. Hoekstra, T. D. Kitching, Y. Mellier, 

L. Miller, J. Coupon, et al., MNRAS 433, 3373 (2013), 
1303.1806. 

[34] The Dark Energy Survey Collaboration, ArXiv Astrophysics e- 
prints (2005), astro-ph/0510346. 

[35] B. Elaugher, International Journal of Modem Physics A 20 , 
3121 (2005). 

[36] H. T. Diehl and for Dark Energy Survey Collaboration, Physics 
Procedia37, 1332 (2012). 

[37] B. L. Elaugher, T. M. C. Abbott, R. Angstadt, J. Annis, M. L. 
Antonik, J. Bailey, O. Ballester, J. P Bernstein, R. A. Bernstein, 

M. Bonati, et al. (2012), vol. 8446 ofProc. SPIE, p. 11. 

[38] B. Elaugher, H. T. Diehl, K. Honscheid, T. M. C. Abbott, 
O. Alvarez, R. Angstadt, J. T. Annis, M. Antonik, O. Ballester, 

L. Beaufore, et al., ArXiv e-prints (2015), 1504.02900. 

[39] C. Chang, V. Vikram, B. Jain, D. Bacon, A. Amara, M. R. 
Becker, G. Bernstein, C. Bonnett, S. Bridle, D. Brout, et al., 
LE15153 (2015), 1505.01871. 

[40] A. Joyce, B. Jain, J. Khoury, and M. Trodden, Physics Reports 
568 , 1 (2015), 1407.0059. 

[41] N. Kaiser and G. Squires, ApJ 404 , 441 (1993). 

[42] M. Bartelmann, R. Narayan, S. Seitz, and P. Schneider, ApJ 
464 , LI 15 (1996), astro-ph/9601011. 

[43] S. L. Bridle, M. P Hobson, A. N. Lasenby, and R. Saunders, 
MNRAS 299 , 895 (1998), astro-ph/9802159. 

[44] J. Liu, A. Petri, Z. Haiman, L. Hui, J. M. Kratochvil, and 

M. May, Phys.Rev.D 91, 063507 (2015), 1412.0757. 

[45] L. Van Waerbeke, MNRAS 313, 524 (2000), astro-ph/9909160. 

[46] P. Simon, M. Hetterscheidt, M. Schirmer, T. Erben, P Schnei¬ 
der, C. Wolf, and K. Meisenheimer, A&A 461, 861 (2007), 
astro-ph/0606622. 

[47] C. Ngeow, J. J. Mohr, T. Alam, W. A. Barkhouse, C. Beldica, 
D. Cai, G. Danes, R. Plante, J. Annis, H. Lin, et al. (2006), vol. 
6270 of Proc. SPIE, p. 23, astro-ph/0608246. 

[48] 1. Sevilla, R. Armstrong, E. Bertin, A. Carlson, G. Danes, S. De- 


sai, M. Gower, R. Gmendl, W. Hanlon, M. Jarvis, et al., ArXiv 
e-prints (2011), 1109.6741. 

[49] S. Desai, R. Armstrong, J. J. Mohr, D. R. Semler, J. Liu, 
E. Bertin, S. S. Allam, W. A. Barkhouse, G. Bazin, E. J. 
Buckley-Geer, et al., ApJ 757, 83 (2012), 1204.1210. 

[50] J. J. Mohr, R. Armstrong, E. Bertin, G. Danes, S. Desai, 
M. Gower, R. Gruendl, W. Hanlon, N. Kuropatkin, H. Lin, et al. 
(2012), vol. 8451 of Proc. SPIE, 1207.3189. 

[51] E. Bertin, in Astronomical Data Analysis Software and Systems 
XV, edited by C. Gabriel, C. Arviset, D. Ponz, and S. Enrique 
(2006), vol. 351 of Astronomical Society of the Pacific Confer¬ 
ence Series, p. 112. 

[52] E. Bertin, Y. Mellier, M. Radovich, G. Missonnier, P. Dide- 
lon, and B. Morin, in Astronomical Data Analysis Software and 
Systems XI, edited by D. A. Bohlender, D. Durand, and T. H. 
Handley (2002), vol. 281 of Astronomical Society of the Pacific 
Conference Series, p. 228. 

[53] E. Bertin and S. Arnouts, A&AS 117 , 393 (1996). 

[54] E. Bertin, in Astronomical Data Analysis Software and Systems 
XX, edited by 1. N. Evans, A. Accomazzi, D. J. Mink, and A. H. 
Rots (2011), vol. 442 of Astronomical Society of the Pacific 
Conference Series, p. 435. 

[55] N. Benitez, ApJ 536 , 571 (2000), astro-ph/9811189. 

[56] D. Coe, N. Benitez, S. E. Sanchez, M. Jee, R. Bouwens, and 
H. Eord, AJ 132, 926 (2006), astro-ph/0605262. 

[57] C. Sanchez, M. Carrasco Kind, H. Lin, R. Miquel, E. B. 
Abdalla, A. Amara, M. Banerji, C. Bonnett, R. Brunner, 

D. Capozzi, et al., MNRAS 445 , 1482 (2014), 1406.4407. 

[58] G. D. Coleman, C.-C. Wu, and D. W. Weedman, ApJS 43, 393 
(1980). 

[59] A. L. Kinney, D. Calzetti, R. C. Bohlin, K. McQuade, 
T. Storchi-Bergmann, and H. R. Schmitt, ApJ 467, 38 (1996). 

[60] G. Bmzual and S. Chariot, MNRAS 344, 1000 (2003), astro- 
ph/0309134. 

[61] C. Bonnett, ArXiv e-prints (2013), 1312.1287. 

[62] P. Graff and E. Eeroz, SkyNet: Neural network training tool 
for machine learning in astronomy. Astrophysics Source Code 
Library (2013), 1312.007. 

[63] E. S. Sheldon, MNRAS 444 , L25 (2014), 1403.7669; down¬ 
loaded at https : //github. com/esheldon/ngmix. 

[64] G. de Vaucouleurs, Annales d’Astrophysique 11 , 247 (1948). 

[65] D. W. Hogg and D. Lang, PASP 125, 719 (2013), 1210.6563. 

[66] A. P. Dempster, N. M. Laird, and D. B. Rubin, JOURNAL 
OP THE ROYAL STATISTICAL SOCIETY, SERIES B 39, 1 
(1977). 

[67] J. Zuntz, T. Kacprzak, L. Voigt, M. Hirsch, B. Rowe, and S. Bri¬ 
dle, MNRAS 434 , 1604 (2013), 1302.0183; downloaded at 
https://bitbucket.org/j oezuntz/imSshape/. 

[68] A. Refregier, T. Kacprzak, A. Amara, S. Bridle, and B. Rowe, 
MNRAS 425, 1951 (2012), 1203.5050. 

[69] T. Kacprzak, J. Zuntz, B. Rowe, S. Bridle, A. Refregier, 
A. Amara, L. Voigt, and M. Hirsch, MNRAS 427, 2711 (2012), 
1203.5049. 

[70] P. Melchior, E. Suchyta, E. Huff, M. Hirsch, T. Kacprzak, 

E. Rykoff, D. Gruen, R. Armstrong, D. Bacon, K. Bechtol, 
et al., MNRAS 449 , 2219 (2015), 1405.4285. 

[71] E. S. Rykoff, E. Rozo, M. T. Busha, C. E. Cunha, 

A. Pinoguenov, A. Evrard, J. Hao, B. P. Koester, A. Leauthaud, 

B. Nord, et al., ApJ 785 , 104 (2014), 1303.3562. 

[72] M. T. Busha, R. H. Wechsler, M. R. Becker, B. Erickson, 
and A. E. Evrard, in American Astronomical Society Meeting 
Abstracts (2013), vol. 221 of American Astronomical Society 
Meeting Abstracts, p. 341.07. 

[73] V. Springel, MNRAS 364, 1105 (2005), astro-ph/0505010. 



20 


[74] M. Crocce, S. Pueblas, and R. Scoccimarro, MNRAS 373 , 369 
(2006), astro-ph/0606505. 

[75] A. Lewis and S. Bridle, Phys.Rev.D 66, 103511 (2002), astro- 
ph/0205436. 

[76] C. Conroy, R. H. Wechsler, and A. V. Kravtsov, ApJ 647, 201 
(2006), astro-ph/0512234. 

[77] R. M. Reddick, R. H. Wechsler, J. L. Tinker, and P S. Behroozi, 
ApJ 771, 30 (2013), 1207.2160. 

[78] H. Aihara, C. Allende Prieto, D. An, S. F. Anderson, 
E. Aubourg, E. Balbinot, T. C. Beers, A. A. Berlind, S. J. Bick- 
erton, D. Bizyaev, et al., ApJS 193 , 29 (2011), 1101.1559. 

[79] J. A. Newman, M. C. Cooper, M. Davis, S. M. Eaber, A. L. 
Coil, P. Guhathakurta, D. C. Koo, A. C. Phillips, C. Conroy, 
A. A. Dutton, et al., ApJS 208 , 5 (2013), 1203.3192. 

[80] R. M. Szepietowski, D. J. Bacon, J. P. Dietrich, M. Busha, 

R. Wechsler, and P. Melchior, MNRAS 440 , 2191 (2014), 
1306.5324. 

[81] M. R. Becker, MNRAS 435, 115 (2013). 

[82] Y. Utsumi, S. Miyazaki, M. J. Geller, I. P. Dell’Antonio, 

M. Oguri, M. J. Kurtz, T. Hamana, and D. G. Eabricant, ApJ 
786 , 93 (2014), 1304.4656. 

[83] E. S. Rykoff, B. P. Koester, E. Rozo, J. Annis, A. E. Evrard, 

S. M. Hansen, J. Hao, D. E. Johnston, T. A. McKay, and R. H. 
Wechsler, ApJ 746 , 178 (2012), 1104.2089. 

[84] M. Cautun, R. van de Weygaert, B. J. T. Jones, and C. S. Erenk, 
MNRAS 441 , 2923 (2014), 1401.7866. 

[85] I. Szapudi, A. Kovacs, B. R. Granett, Z. Erei, J. Silk, W. Burgett, 
S. Cole, P. W. Draper, D. J. Earrow, N. Kaiser, et al., MNRAS 
450 , 288 (2015), 1405.1566. 

[86] H. T. Diehl, T. M. C. Abbott, J. Annis, R. Armstrong, L. Baruah, 
A. Bermeo, G. Bernstein, E. Beynon, C. Bruderer, E. J. 
Buckley-Geer, et al. (2014), vol. 9149 of Proc. SPIE, p. 0. 

[87] A. J. Ross, W. J. Percival, A. G. Sanchez, L. Samushia, S. Ho, 
E. Kazin, M. Manera, B. Reid, M. White, R. Tojeiro, et al., 
MNRAS 424 , 564 (2012), 1203.6499. 

[88] S. Ho, A. Cuesta, H.-J. Seo, R. de Putter, A. J. Ross, M. White, 

N. Padmanabhan, S. Saito, D. J. Schlegel, E. Schlafly, et al., 
ApJ 761, 14 (2012), 1201.2137. 

[89] http://WWW.darkenergysurvey.org 

[90] We adopt the “M0DEST_CLASS” classifier, which is a new clas¬ 
sifier used for SVAl Gold that has been developed empirically 
and tested on DES imaging of COSMOS fields with Hubble 
Space Telescope ACS imaging. 

[91] The open source code can be downloaded at: https:// 
github.com/esheldon/ngmix 

[92] The open source code can be downloaded at: https:// 
bitbucket.org/j oezuntz/imSshape/ 


Appendix A: Foreground sample selection 

As discussed in Sec. IIIB, we consider two factors that can 
affect the selection of our foreground sample - spatial varia¬ 
tion in depth and spatial variation in seeing. If not taken care 
of, these effects will result in apparent spatial variation of the 
foreground galaxy number density that is not due to the cos¬ 
mological clustering of galaxies. Below we describe tests for 
each of these and determine a set of selection criteria based on 
the analysis. 


1. Depth variation 

Spatial variation in the depth of the images can cause the 
apparent galaxy number density to vary, as more or less galax¬ 
ies survive the detection threshold. We would like to con¬ 
struct a foreground galaxy sample which minimizes this vary¬ 
ing depth effect. A simple solution is to place a magnitude 
selection slightly shallower than the limiting magnitude in all 
of the areas considered, so that the sample is close to complete 
in that magnitude range. 

We find that in our area of interest, with a magnitude selec¬ 
tion at i <22, we have 97.5% of the area that is complete to 
this magnitude limit. We use the lOcT galaxy limiting magni¬ 
tude to define depth, which is a conservative measure for the 
completeness, as we detect many more galaxies below lOcT. 
The detail methodology of estimating the limiting magnitude 
of the data is described in Rykoff et al. (in preparation). The 
2.5% slightly shallower is not expected to yield significant 
change in our results. 


2. Seeing variation 

Spatial variation in seeing can lead to spatial variation in 
apparent galaxy number density, as large seeing leads to less 
effective star-galaxy separation as well as higher probability 
of blending in crowded fields. To test this, we first select a 
foreground sample with i < 22 and 0.1 < z < 0.5 according to 
Sec. IIIB. Then we look at the correlation between the galaxy 
number density in this foreground sample and the average see¬ 
ing values at these locations, both calculated on a grid of 5 x5 
arcmin^ pixels without smoothing. Fig. (13) shows the galaxy 
number density versus seeing. The black data points show the 
mean and standard deviation (multiplied by 10 for easy visu¬ 
alisation) of the scatter plot in 15 seeing bins. There is a small 
anti-correlation between these two values at the 6% level. This 
is at an acceptable level for us to continue the analysis without 
masking out the extreme high/low seeing regions. 


Appendix B: Correcting for systematic contamination using 
PCA 

As shown in Sec. VII, we can use Eq. (18) to check for any 
outstanding systematic contamination in our Ke map and its 
correlation with the Kg map. Here we present a general treat¬ 
ment to correct for these systematic contaminations, similar 
to that used in Ross et al. [87] and Ho et al. [88]. 

Assume that our measured Ke map is a linear combination 
of the true and some small coefficient times the 

systematics maps {Mi} that can potentially contaminate the 
Ke maps (e.g.seeing, PSF ellipticity). That is 

N 

— NE'^true (^ 1 ) 

i 

where we have a total of N systematics maps. Similarly, we 
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FIG. 13: Galaxy number density as a function of the seeing in the 
area of consideration. The black line shows the mean and standard 
deviation (multiplied by 10 for easy visualisation) of the scatter plot 
in 15 seeing bins. 
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FIG. 14: The systematics map for Ke (left) and Kg (right) shown 
is compiled using a linear combination of 20 principal components 
extracted from the systematics maps listed in Table II. 


have the expression for the measured Kg in the same way 


N 

— ^g,true PiMi , 

i 

(B2) 

where ft is the linear coefficient in this case. 

Assuming the true maps are uncorrelated with the system¬ 
atics maps, we have 

(^£',true^i) — 

(B3) 

(^g,trueMi) = 0. 

(B4) 

Correlating the measured Ke with a single systematics map 
gives 

N 

{keMj) = {(£aiMi)Mj). 
i 

(B5) 

We can construct a set of systematics maps that are uncorre¬ 
lated between each other, or (MiMj^i) = 0, and then extract 
all the coefficients at from the observables as follows: 

{keMj) = ajiMjMj); 


„ i^EMj) , 

' {MjMjY 


KEtn>c = KE f j M,. 

i.true £ 

(B6) 

And similarly for Kg, we have 



(B7) 


To construct a set of systematics maps {Mi} uncorrelated be¬ 
tween each other from a set of systematics maps correlated 
with each other {M-} (i.e. those listed in Table II), we in¬ 
voke the Principal Component Analysis (PCA) method. In 
this case, each of the pixelated maps, after normalizing by its 
scatter, {M-} form a data vector, and the extracted eigenvec¬ 
tors form a orthogonal basis set, which we can use as {Mi}. 
We find that the principal component maps correspond strik¬ 
ingly to physical properties of the data. Fig. (14) shows the 
systematics maps corresponding to Ke and main sample Kg 
extracted using this PCA method, or the second terms on the 
right-hand-side of Eq. (B6) and Eq. (B7). We find that the 
main contributions come from large-scale structures and are at 
a very low level compared to the original maps (see Eig. (2)). 
We subtract these systematics maps from the original Ke and 
Kg maps according to Eq. (B6) and Eq. (B7). The Pearson 
correlation coefficient changes by 3.5% relative to the origi¬ 
nal pKEKg measured in Sec. VI, suggesting the contamination 
to the cross-correlation coefficient is not significant. 

It is worth noting that there are a few assumptions that go 
into the calculation above, which need to be accounted for 
when interpreting these results. Eirst, we have assumed that 
the systematic maps have no correlation with the true Ke and 
Kg maps. Eor a large enough area, this should be true, but for 
small maps we can expect some correlation just by chance. 
Hence the quantitative “improvement” we get in the Pearson 
correlation coefficient must be carefully checked with simu¬ 
lations with larger area than used here. Second, since the 
method is based on PCA, the effectiveness of the correction 
depends on finding the important systematics maps that can 
contribute linearly to the contamination. That is, if the sys¬ 
tematics come from a non-linear combination of the various 
maps (e.g. multiplication of two maps), then one would not 
automatically correct for it without putting in this correct non¬ 
linear combination of maps in the first place. 
















